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Non-compact QED3 with four- component fermion flavor content Nf > 2 is studied numerically near the chiral 
limit to understand its chiral symmetry breaking features. We monitor discretization and finite size effects on 
the chiral condensate by simulating the model at different values of the gauge coupling on lattices ranging in size 
from 10^ to 50^ Our upper bound for the dimensionless condensate (3^{'i>'i>) in the Nf = 2 case is 5 x 10 ^. 
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1. Introduction 

Over the last few years, QED3 has attracted a 
lot of attention because of its potential applica- 
tions to models of high T^. superconductivity. It is 
believed to be confining and exhibit features such 
as dynamical mass generation when the number 
of fermion flavors Nf is smaller than a critical 
value Nfc- It is therefore an interesting and chal- 
lenging model and an ideal laboratory in which 
to study more complicated gauge field theories. 

We are considering the four-component for- 
mulation of QED3 where the Dirac algebra is 
represented by the 4x4 matrices 70, 71 and 
72. This formulation preserves parity and gives 
each spinor a global U{2) symmetry generated 
by 1, 73,75 and 17375; the full symmetry is then 
U{2Nf). If the fermions acquire dynamical mass 
the U(2Nf) symmetry is broken spontaneously to 
U{Nf) X U{Nf) and 2iV| Goldstone bosons ap- 
pear in the particle spectrum. 

At the present time the issue of spontaneous 
chiral symmetry breaking in QED3 is not very 
well understood. Studies based on Schwinger- 
Dyson equations (SDEs) using the photon prop- 
agator derived from the leading order 1/Nf ex- 
pansion suggested that for Nf less than some 
critical value Nfc the answer is positive with 
Nfc ~ 3.2 0. Apparently, for Nf > Nfc, the at- 
tactive interaction between a fermion and an an- 
tifermion due to photon exchange is overwhelmed 



by the fermion screening of the theory's electric 
charge. More recent studies which treat the ver- 
tex consistently in both numerator and denom- 
inator of the SDEs have found a value for Nfc 
either in agreement with the original study or 
slightly higher, with Nfc ~ 4.3 |||]. Finally an 
argument based on a thermodynamic inequality 
has predicted Nfc < | ||]. There have also been 
numerical attempts to resolve the issue via lat- 
tice simulations. The obvious advantage in this 
case is that one can study any Nf without any 
assumption concerning the convergence of expan- 
sion methods. However, the principal obstruction 
to a definitive answer has been large finite volume 
effects resulting from the presence of a massless 
photon. Numerical studies of the quenched case 
have shown that chiral symmetry is broken [Q, 
whereas in the case of simulations with dynam- 
ical fermions opinions have divided on whether 
Nfc is finite and « 3 or whether chiral symme- 
try is broken for all Nf Q . Recent studies of the 
Nf = 1 model on small lattices appear in . 

In our study we used the Hybrid Monte Carlo 
algorithm to simulate the non-compact version 
of lattice QED3 with staggered fermions. In the 
continuum it corresponds to the four-component 
spinor formulation of the model 0. We imple- 
ment even-odd partitioning which implies that a 
single flavor of one-component staggered fermions 
can be simulated, which corresponds to Nf = 2 
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in the continuum limit. 

2. Results 

In this section wc discuss the resuhs of lattice 
simulations of QED3 with Nf > 2. In our study 
we tried to detect and control the various draw- 
backs of the lattice method: (i) The lattice itself 
distorts continuum space-time physics consider- 
ably unless the lattice spacing a can be chosen 
small compared to the relevant physical wave- 
lengths in the system, (ii) the size of the lat- 
tice must be large compared to the dynam- 
ically generated correlations in the system; and 
(iii) the chiral limit can only be studied by simu- 
lating light fermions of mass mo in lattice units. 
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Figure 2. Dimensionless condensate vs. dimen- 
sionless bare mass for N i = 2 at different values 
of the coupling /? and constant physical volume 
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Figure 1. Dimensionless condensate vs. 
dimensionless bare mass /3mo for Nf = 2, 4, 8, 16, 
(3 — 0.6 on a 16^ lattice. 

In Fighl we plot the dimensionless chiral con- 
densate p^{^'^) vs. the dimensionless bare mass 
(3mo (where P = for Nf = 2,4,8,16. The 
coupling (3 = 0.6 and the lattice volume is 16'^. 
As Nf increases the chiral condensate decreases 
(for Too > 0) because the interaction between the 
fermion and the antifermion is screened. How- 
ever, as mo — > all the curves tend to pass 
smoothly through the origin. This motivated us 
to study in more detail the pattern of chiral sym- 
metry breaking at small Nf on larger volumes 
near the chiral limit. 

In order to check whether our lattice data are 
characteristic of the continuum limit we plot in 



FigJI vs. /3too for Nf = 2 and couphng 

(3 = 0.45,0.60,0.75,0.90. To disentangle the lat- 
tice discretization effects from the finite size ef- 
fects we keep the volume in physical units (L//3)^ 
constant. It can be inferred from the graph that 
discretization effects are small for /? > 0.60, be- 
cause the data almost fall on the same line within 
the resolution of our analysis. 

In Fig.^ we present our results for the chiral 
condensate vs. bare mass for Nf = 2 and P = 0.6 
on lattice sizes varying from 8'^ to 48"^. We infer 
that finite size effects become small for L > 24 
and all the lines tend to pass smoothly through 
the origin. Our analysis of meson masses and 
susceptibilities in scalar and pseudoscalar chan- 
nels showed that these quantities suffer from very 
strong finite size effects and therefore did not al- 
low us to reach such a definitive conclusion [^. 

Next we discuss the results from Nf — 2 sim- 
ulations at (3 — 0.75. The (3 = 0.75 data set is 
closer to the continuum limit than the data ex- 
tracted at /3 = 0.60. However, particular care 
is required because weak coupling data are very 
sensitive to finite size effects and accurate mea- 
surements require simulations on large lattices. 
In Fig.|| we present the results for the chiral con- 
densate vs. fermion bare mass extracted from 
simulations with lattice sizes ranging from 10"^ to 
50^^. These simulations were performed very close 
to the chiral limit, i.e. with mo < 0.005. We can 



3 



see from the figure the finite size effects are under 
relatively good control and the data tend to pass 
smoothly through the origin. Therefore, we con- 
clude that for Nf = 2 < 5 x 10-^ which 
is a strong indication that QED3 may be chirally 
symmetric for Nf > 2. 
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g ^ 0, on different volumes in order to detect 
and control finite size effects and near the chiral 
limit niQ — > 0. 

Our upper bound for the condensate in the 
Nf = 2 case is < 5 x 10"^ and all the 

hues of /3^(^\E') vs. (3nio tend to pass smoothly 
through the origin which may imply that chiral 
symmetry is restored for Nf > 2. We are con- 
tinuing this study to check if chiral symmetry is 
broken in the case of Nf — 1. 
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Figure 3. Dimensionless condensate vs. dimen- 
sionless bare mass for Nf = 2, f3 = 0.6 and lattice 
sizes 83,163,243,323,483. 
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Figure 4. Condensate vs. bare mass for Nf — 2, 
(3 = 0.75 and lattice sizes 103, 203, 303, 403 and 
503. 



3. Conclusions 

In our study of QED3 with Nf > 2 we at- 
tempted to establish whether chiral symmetry is 
broken or not by studying the behavior of the 
chiral condensate close to the continuum limit 
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